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We study the two-body decays of B mesons to K and a charmonium state, X C c, in a sample of 
210.5 fb _1 of data from the BABAR experiment. We perform measurements of absolute branching 
fractions B(B ± -> K^X cB ) using a missing mass technique, and report several new or improved 
results. In particular, the upper limit B(B ± -> K ± X(3872)) < 3.2 x 10" 4 at 90% CL and the 
inferred lower limit B(X (3872) — + J/ipir + -K~) > 4.2% will help in understanding the nature of the 
recently discovered A" (3872). 

PACS numbers: 13.25.Hw, 14.40.Gx 



Several exclusive decays of B mesons of the form 
B^ — » K^X cS (where X cS is one of the charmonium 
states rj c , J/ip , \c0 , Xci , rj'c > ; V'" )■> have been observed 
by reconstructing the charmonium state from its decay 
to some known final state, / In principle, such B 

decays provide a direct probe of charmonium properties 
since the phase space is large for all known states and 
all should be produced roughly equally, in the absence 
of a strong selection rule [3j. However with this tech- 
nique only the product of the two branching fractions 
B(B ± — > K ± X cS )xB(X cE — > /) is measured, thereby 
reducing the precision of £>(B ± — » K ± X cS ) when the 
daughter branching fraction is poorly known. 

We describe here a complementary approach, based on 
the measurement of the kaon momentum spectrum in the 
B center-of-mass frame, where two-body decays can be 
identified by their characteristic monochromatic line, al- 
lowing an absolute determination of B(B ± — > K^X cE ). 
Knowledge of the B center-of-mass system is obtained 
by exclusive reconstruction of the other B meson from 
a T(AS) decay. In addition to obtaining new informa- 
tion on known charmonium states, this method is used 
to search for the A (3872) state, recently observed in 
B ± -> # ± X(3872) decays by Belle and BABAR 0, 
in the subsequent decay A(3872) — > J j 't/m + ir~ . The 
same method allows a search for charged partners of the 
X(3872) in B° decays, independent of the X(3872) ± de- 
cay mode. The nature of AT(3872) resonance is still un- 
clear, different interpretations 6] have been proposed but 
more experimental data will be needed to discriminate 
between them. 

For this analysis we use a data sample of 210.5 fb _1 
integrated luminosity, corresponding to 231.8 x 10 6 BB 
pairs. The data have been collected with the BABAR de- 
tector at the SLAC PEP-II asymmetric-energy collider, 
where 9 GeV electrons and 3.1 GeV positrons collide at 
a center-of-mass energy 10.58 GeV, corresponding to the 
mass of the T(4S) resonance. A detailed description of 
the BABAR detector can be found in ■ Charged tracks 
are reconstructed with a 5 layer silicon vertex tracker 
(SVT) and a 40 layer drift chamber (DCH), located in 
a 1.5 T magnetic field generated by a superconducting 



solenoid. The energy of photons and electrons is mea- 
sured with an electromagnetic calorimeter made up of 
CsI(Tl) crystals. Charged hadron identification is done 
with ionization measurements in the SVT and DCH and 
with an internally reflecting ring imaging Cherenkov de- 
tector. The instrumented flux return of the solenoid is 
used to identify muons. 

The analysis is performed on a sample of events where 
a B meson is fully reconstructed (B recon ). For these 
events, the momentum of the other B (B S i gna i) can be 
calculated from the momentum of B recon and the beam 
parameters. We select events with a K^ not used for 
the reconstruction of B recon and calculate its momentum 
(Pk) in the B S i gna i center of mass system. 

B recon mesons are reconstructed in their decays to ex- 
clusive D(*>H final states, where H is one of several com- 
binations of ir^, K^, 7T° and Kg hadrons; a detailed de- 
scription of the method can be found in [8j . 

The number of B^ events in the data is determined 
with a fit to the distribution of the beam energy substi- 
tuted mass niEs = ^B^, M /4 — p 2 B , where Eqm is the 
total center-of-mass energy, determined from the beam 
parameters, andps is the measured momentum of B recon 
in the center-of-mass frame. The fit function is the sum 
of a Crystal Ball function describing the signal and 
an ARGUS function ^} f° r each background component 
(e + e~ — > qq where q is u, d, s or c or misreconstructed 
Bs), the relative weights of which are obtained from a 
Monte Carlo simulation (MC), while the total normal- 
ization factor is determined from the data. A total of 
378580 ±1110 events with a fully reconstructed B^ 1 is 
obtained. 

Fifteen variables related to the B recon decay character- 
istics, its production kinematics, the topology of the full 
event, and the angular correlation between B recon and 
the rest of the event are used in a neural network (NN1) 
to reduce the large background, mainly due to non-B 
events. The network has 80% signal efficiency while re- 
jecting 90% of the background. The tties distribution 
after this selection is shown in Fig. ^ Only events with 
5.275 < rriEs < 5.285 GeV/c 2 are used in the analysis. 

We now consider only tracks not associated with 
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FIG. 1: The rriES distribution of all B recon after the NN1 
selection. The solid line represents a fit described in the text; 
cumulative background contributions from e + e~ — > qq where 
q is u, d, s or c (dashed line), _B° (dotted line), B (dash-dot 
line) events are shown. The arrows indicate the cuts used in 
the analysis (see text). 



B recon . Most K ± produced in B^ decays originate from 
D mesons and their spectrum, although broad, peaks at 
low pk- In the B^ rest frame, these if ± are embedded 
in a "minijet" of D decay products, while signal re- 
coil against a massive (3-4 GeV/c 2 ) state and therefore 
tend to be more isolated. A second neural network (NN2) 
rejects background from secondary K , by using fifteen 
input variables describing the energy and track multi- 
plicities measured in the if hemisphere, the sphericity 
of the recoil system, and the angular correlations be- 
tween the K ± and the recoil system. These variables 
have been chosen to be independent of the particular de- 
cay topology of the recoil system. Since the topology of 
the event changes with the recoil mass, we have consid- 
ered separately two recoil mass regions in the training 
of this neural network: the "high-mass" region, corre- 
sponding to 1.0< pk <1.5 GeV/c and the "low-mass" 
region, for 1.5< pk <2.0 GeV/c. The signal training 
sample is S ± — » K ± X c5 MC simulation while the back- 
ground sample consists of simulated K from D meson 
decays in the same momentum range. The chosen cuts 
on the NN2 outputs correspond to 85% signal efficiency; 
the background rejection factor varies between 2.5 in the 
A(3872) and tp' region and 1.5 in the J/ip region. The 
selection criteria are optimized for MC signal significance 
with the high-mass region blinded. 

The kaon momentum distribution shows a series of 
peaks due to the two-body decays B ± — > K ± X c5 cor- 
responding to the different X c£ masses, superimposed 
on a smooth spectrum due to if ^ coming from multi- 
body B ± decays, or non-S ± background. The mass 
of the X c5 state (mi) can be calculated directly from 
Pk using mx = y rn% + m 2 K — 2EKrriB, where raj and 
rriK are the B^ and if ± masses and Ek is the if ± en- 
ergy. The resonance width Tx can be obtained from 
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FIG. 2: Kaon momentum spectrum for the (a) low-mass and 
(b) high-mass regions. The lines represent the fit described in 
the text. Arrows show the the expected positions of known 
charmonium states. 



the Breit-Wigner width of the peak in the pjcspectrum 
r^-, obtained after deconvolution with the momentum 
resolution function, using Tx — TkPkitib/itix, where 
Pk =pk/E K - 

We determine the number of S ± — » K ± X cS events 
(Nx) from a fit to the pk distribution. The branching 
fraction for the decay channel is calculated as: 

B(B ± -» ^X, " A X 



ex ■ N B ■ 



where tx is the efficiency determined from the MC and 
Nb the number of B mesons in the sample. An alter- 
native method, which we use to improve the branching 
fraction measurement in the case of r\ c , is to normalize 
to the channel B^ — * K J/tp, which is well- measured in 
the literature [TJ, according to: 



B{B d 



K ± X r 



N x 



N 



7/i/> e X 



B(B ± -► K^J/ijj). 



In this relative measurement, the systematic errors that 
are common to both resonances cancel in the ratio. The 
two methods are combined to extract B(B ± — * K^rjc), 
taking into account the correlations between them. 

We fit the pk spectrum using an unbinned maximum 
likelihood method. The background is well modeled by a 
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TABLE I: Event yields and absolute branching fractions 
B(B ± — > K X c s) from the fits to the Pk spectrum. The 
first error is statistical, the second systematic and B upper 
limits are given at 90% CL, taking into account the 9.0% sys- 
tematic error. The last column shows the signal statistical 
significance a, derived from the fit likelihood assuming sig- 
nal events L(0): a = yj —2 log L(0). For the r\ c , both results 
for absolute and relative measurement, and their combination, 
are reported (see text). 



Particle 


Yield 


S(HT 4 ) 


a 


Vc 


273±43 


8.4±f.3±0.8 


7.3 


r/ c relative 




f0.6±2.3±0.4±0.4 




rj c combined 




8.7±1.5 




J/i> 


259±41 


8.1±f.3±0.7 


6.9 


XcO 


9±21 


<1.8 




Xcl 


227±40 


8.0±f.4±0.7 


6.0 


Xc2 


0±36 


<2.0 




Vc 


98±52 


3.4±f.8±0.3 


1.8 


V>' 


139±44 


4.9±1.6±0.4 


3.2 


r 


99±69 


3.5±2.5±0.3 


1.4 


A(3872) 


f5±39 


< 3.2 





third degree polynomial and each signal is a Breit-Wigner 
function folded with a resolution function. The masses 
and widths of the r/ c and r)' c mesons are left free; all others 
are fixed to values from reference The resolution 

function has two parts: a Gaussian with a varying from 
6 MeV/c at p K ~ 1.1 GeV/c to 12 MeV/c at p K ^ 1.7 
GeV/c describes the 72.5% of the signal where B recon 
is correctly reconstructed; if B recon is incorrect, but has 
mES within our range, the px resolution is a bifurcated 
Gaussian with a = 78 and 52 MeV/c on the left and right 
hand side of the peak respectively. 

The spectrum in the low-mass region is expected to 
exhibit two peaks, at px — 1.683 GeV/c corresponding 
to the J/ip , and at px = 1.754 GeV/c for the t] c meson. 
These two peaks are clearly seen in Fig. Ufa); both have 
a significance of ~ 7a. The number of events under each 
peak obtained from the fit is N(J/ip) = 259 ± 41 and 
N(i lc ) = 273 ± 43. 

The spectrum in the high-mass region is fitted with 
a background and seven signal functions, corresponding 
to the following states: ip' , x c0 , Xci, Xc2, ip", v' c and 
X(3872). The resulting fit is shown in Fig. HTb), with 
the yields given in Tabled The h c charmonium state lies 
near the Xci, and it is difficult to distinguish the peaks 
from these two decays. A fit including the h c yields a 
number of h c events consistent with zero, and a fit per- 
formed with free Xcl mass and width gives values con- 
sistent with a narrow Xci , therefore we have no evidence 
for h c production. 

Several sources of systematic error affecting these mea- 
surements have been evaluated. The relative errors on 
absolute measurements are the same for all states; many 
of these cancel partially in relative measurements, and all 
are summarized in Table ITT1 "B counting" refers to un- 



TABLE II: Summary of systematic errors in percent for ab- 
solute and the J/ip : 7/ c relative measurement. 



Source 


a i i j_ /fry \ 

Absolute (%) 


J/ip -.r/c (%) 


B counting 


4.5 





Mass scale 


1 


1 


Background model 


3.5 


1.7 


Resolution model 


2.3 


1.0 


K ± reconstruction 


1.3 





K identification 


5 


1 


B mass selection 


0.5 





NN1 selection 


2.2 


2.0 


NN2 selection 


3.2 


1.0 


Total 


9.0 


3.3 



certainties in the fit parametrization used to determine 
the number of fully reconstructed B^ econ . It is one of 
the largest errors in absolute measurements, and cancels 
in ratios. The mass scale is verified to a precision of 1.5 
MeV/c in pk by floating the masses of the well- measured 
J/ip i Xci and ip' peaks; we assign a systematic error 
corresponding to this shift. We also consider variations 
in the background and signal model parametrizations, 
which partially cancel in the case of ratios. Errors in 
the track reconstruction and identification efficiency 
are evaluated by comparing data and MC control sam- 
ples. The systematic error in the NN1 and NN2 selections 
is evaluated by comparing efficiencies and distributions 
in data and MC, and studying efficiency variation with 
Pk- We verified that the NN2 selection is not depen- 
dent on visible energy or multiplicity of the recoil part 
of the B meson decay. Adding in quadrature, the total 
relative error on an absolute measurement is 9.0%. The 
total is reduced to 3.3% for the relative measurement of 
J/ip and r/c , and to 5.9% for states in the high-mass 
region relative to J/ip . For the extraction of relative 
branching fractions, an additional 4% error, labeled (ext) 
in the following, comes from the present uncertainty of 
B(B ± -> K ± J/i[>) = (10.0 ± 0.4) x 10" 4 [13. 

In the high-mass region, clear signals are found for 
Xci and ip' (with significance 6.0 and 3.2cr respectively), 
an excess of events is present for rj' c and ip" [12j , while no 
signal is found for XcO , Xc2 and AT(3872). The branching 
fractions and upper limits are summarized in Table [I] 

In the low-mass region, our J/ip measurement is consis- 
tent with the world average. From the r\ c and J/ip yields 
and the reference branching fraction we can derive the 
result with the relative measurement method B(B ± —> 
K^)^ = (10.6±2.3(stat)±0.4(sys)±0.4(ext))xl0~ 4 . 
We combine this result with the absolute measurement 
of Table |l] taking the correlated errors into account, to 
obtain B{B ± -> K^rjc) = (8.7 ± 1.5) x 10~ 4 . 

We obtain from our fits the rj c and rj' c masses and 
widths and find m Vc = 2982 ± 5 MeV/c 2 , r r)c < 43 MeV 
and m n , c = 3639 ± 7 MeV/c 2 , T v>c < 23 MeV, where the 



7 



width limits are both at 90% CL. 

Taking B(B ± -> K ± XJ3872)) < 3.2 x 10~ 4 , and using 
an average of the Belle and BABAR [f| measurements 
of B{B ± -► Ar ± A(3872))xS(X(3872) -> J/^tt+tt") wc 
set a lower limit £(A(3872) -> J/V>7r+7r-)> 4.2% at 90% 
CL. This branching fraction, for which there are not yet 
any predictions, is sensitive to the distribution of charm 
quarks inside the A (3872). A search for charged partners 
of the X(3872) is performed by examining recoiling 
from a sample of 245.6k reconstructed B° decays. No 
signal is seen and we find B{B° K ± X(3872)T) < 
5 x 10~ 4 at 90% CL. 

We combine our B(B ± — > K rj c ) with a previous 
BABAR measurement of B{B ± _-^ K ± r] c )xB(r ]c -> KKn) 
to obtain B(r] c -> KKtt)=(8.5 ± 1.8)%, signif- 
icantly improving the precision of the world average 
Since this branching fraction is used as a refer- 
ence for all r\ c yield measurements, our result will lead 
to more precise rj c partial widths and more stringent 
comparisons with theoretical models. For example, from 
an average of B(Jhb — > "/TjAxB^c — > KKtt) mea- 
sured by Mark-Ill Q, DM2 and BES [II], we ob- 
tain B(J/tp -> 777 c )=(0.79±0.20)%, and usingthe value 
r(?y c -> 77) xB(r? c -> AAvr)=0.48 ± 0.06 keV [ill we cal- 
culate T(rj c — > 77)=(5.6±1.4) keV. Both results are more 
precise than the world average [ill ] . Similarly, we obtain 
B« KKV)=(8±5)% and T(rf c -> 7 7 )=(0.9±0.5)keV. 

In conclusion, a novel technique is used to measure 
directly the absolute branching fractions of the vari- 
ous charmonium states X cg in two-body decays B —> 
K ± X c£ (Table [J. The results for X c£ — rj c , J/ip,ip' are 
in agreement with previous measurements, and the r\ c re- 
sult significantly improves the present world average. Up- 
per limits are set for Xco and \c2 , confirming factoriza- 
tion suppression ^3]- Measurements of B ± — > K ±r q' c 
and B ± — > K ± t/j" branching fractions are reported, al- 
though with poor significance. Upper limits are given for 
A (3872) and for production of a possible charged partner 
in B° decays. 
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